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We show that the coherent mixing of different transverse modes, due to forward scattering of 
carriers by soft impurity- or boundary potentials leads to a nonlinear, asymmetric current response 
of quantum point contacts (QPC). The oscillating contribution to the current is sensitive both to 
driving voltage and to gate voltage in direct analogy to the electrostatic Aharonov-Bohm effect. 

Our calculations are in a good agreement with recent experimental data showing small-scale 
conductivity nonlinearities and asymmetry in QPC. 
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It is well known that the quantization of the electron's transverse momentum in quantum point contacts (QPC) - 
structures, defined in the 2D electron gas of GaAs-AlGaAs interfaces by the electrostatic potential of the gate electrodes 
- is responsible for their nonlinear behaviour both in regard to the gate voltage, V g , and the driving voltage, V s d- The 
Pj^ ■ nonlinearities are due to effective turning on-off the conducting modes inside the contact when eV g or eV s d 0-0] 

exceeds the interlevel spacing of transverse electronic modes in the contact region, AE± ~ Ep/N ~ lmeV (N-number 
of open modes). They were observed experimentally, as differential conductance quantized steps vs. V g or as 
extrema of the conductance vs. V s d derivative (?]]. 

The quantization effects in QPC proved to be stable in regard to the potential variations in the system (e.g., due to 
charged impurities). The reason is that backscattering is needed in order to violate the quantization, while scattering 
on soft potentials present in the system cannot produce a sufficient momentum change |||| ■ This explains the success 
of an adiabatic model of QPC || in spite of significant mode mixture due to forward scattering inside the system ||. 

Nevertheless, these results don't exhaust all the possibilities, and in recent experiments jl0| a more complicated 
picture was observed, which cannot be explained in the framework of the existing approach. 
\ The nonlinearity and asymmetry ( sensitivity to the current direction ) of QPC response was observed at driving 

■ voltages as small as 0.01 mV, both directly in measurements of the differential resistance vs. current and via ac 
current rectification. The effect was clearly seen at 0.3 K and on increasing the temperature till 4.2 K it was smeared 
very similarly to the smearing of the transconductance, dG/dV g . Resistance variation was of order 20 Q, on the 
background of QPC resistance of several kil (QPC conductance ranged from 1 to 10 (2e 2 /ft.)-quanta). 

The purpose of the present paper is to demonstrate that the scattering on soft potentials can play a major role in 

■ the ballistic transport through microconstrictions due to coherent mixing of the transverse modes, their phases being 
q \ defined by V s d and V g , in a direct analogy to the electrostatic Aharonov-Bohm effect. 

The effect exists in a QPC with at least two scatterers (impurities, defects, nonadiabaticities of the contact's shape 
etc). In a system with only soft scatterers a crucial role in its development is played by the mechanism of indirect 
backscattering suggested by Laughton et al. |ll| plays . (It is related to the possibility for a carrier to be switched 
' from the transport mode to non-propagating one even by the forward (small- angle) scattering.) 

We will show, that a qualitative understanding of the measured oscillations and asymmetry of the contact's conduc- 
tance is achieved. The characteristic scales of the oscillations vs. V s d, V g are in good accordance with the experimental 
data as well. 

The model we consider is an adiabatically smooth 2D channel connecting two equilibrium reservoirs, with two 
scatterers inside. The wave function of an electron with energy E in some point of our system can be expanded over 
WKB eigenf unctions for transporting modes || 

r(x,y;E) = '£a«x% l (x,y;E); (l) 



Xm( x >y> E ) = \Jpoc{E)/p m \\{E,x) exp ^i/h ^ dx p m \\ (as'; E)\ 4> m (y; x). (2) 

In this expression 4> rn (y;x) is a transverse eigenfunction with eigenvalue E m ±(x) = (p m ±(x; E)) 2 /2m*; index 
a — 0, 1 numbers the waves propagating from the right (a = 0) or left (a = 1) reservoir. 
The phase gained by a partial wave between the scatterers A,B is 



a n (E) = l/h dx'p n[{ (x';E) = l/h dx'y/2m (E - E n± (x) - e$(x)), (3) 
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and is evidently dependent on the mode number, electron energy and absolute value of the external electrical potential 
in the interscatterer region, Phase differences between different partial waves are thus sensitive to $(x), in a 

direct analogy to the electrostatic Aharonov-Bohm effect, and so is an interference contribution to the current. 
The current through the system is written in a standard way [12] : 



I(V sd ) = e/(irh) J dE [n F (E - m ) - n F (E - fi L )] Tr (f (E)^f (E)) 



(4) 



Here f°(E) is the transfer- matrix for the particle with energy E incident from the right reservoir: ^ ut (E) = 
T° (E^^^E) . N-dimensional vector columns ^f ntOUt contain the expansion coefficients dm-inout Eq.(^) (amplitudes 
of transverse modes in in- and outcoming waves incident from the right or left reservoir respectively) . N is a number 
of conducting modes in the channel. 

The scatterers are described by unitary (2Nx2N) S-matrices [Q. If we introduce the phase gain matrix, Uij = 
Sij exp(i<jj(E)) } then the transfer-matrix is expressed through U and elements of S-matrices of scatterers as follows: 

f°(E) = i%Ut A + P B Uf A Uf B Uf A + . . . = i% (i - Uf A Uf B ) ' Ut A . (5) 

Here t AB ,r AB are (NxN)-matrices which describe the transmission (reflection) to the left (a = 0), or to the right 
(a = I) (see e.g. 13 ) . Summation of higher-order terms, corresponding to multiple passages of the wave through the 
system, allows us to obtain conductance resonances like the ones predicted for QPCs with sharp edges |l4|~|l6|], which 
could contribute to small-scale nonlinearities in the system. This contribution is negigible in our case because of two 
reasons. First, these resonances are due to direct backscattering into the same mode and were never obtained in more 
realistic models with contacts having a smooth shape. Second, the period of the corresponding nonlinearities is much 
smaller than that of oscillations due to intermodc mixing, and they surely would have vanished at the temperature 
of the experiment. 

Keeping thus only the first term of the series, which contains effects of intermode mixing, we find that the current 
is given by ( [|) with 

Tr (T°(E)tf°(E)) =J2H (M) ., . . exp (^(E) - a k (E))) . (6) 

j k 3 J 

The effect we are interested in is contained in phase-sensitive off-diagonal terms of the double sum. 

The exact form of the matrices t cannot be determined without using some specific model of a scatterer, be it 
impurity, contact-bank interface or something else. We are not here concerned with these poorly controllable details, 
and concentrate on the nonlinear conductivity of QPC as a function of V s d and V^f. Nevertheless some general 
remarks are to be made. 

Unitarity condition imposes constraints upon the transition and reflection submatrices. In the general case Jl5| ] 

M) = fi x (i + A^)- 1 ^; (f B %) = 4(1 + Xb)- 1 ^, (7) 

where ux,2 are arbitrary unitary (NxN)-matrices, and A is a diagonal (NxN)-matrix with nonnegative elements. A 
equals to zero if and only if the reflection submatrix is zero, i.e., if the direct backscattering is absent. 

In the latter case both {f^V^j — ^r^jj) — 1> and interference terms are exactly zero, no matter how strong the 

mode mixing is. Thus, as long as number of open modes is not changed (V s d < AVj_), conductance quantization is 
not affected by finite driving voltage at all. This is consistent with numerical calculations by Brataas and Chao ||, 
where mode mixing was shown to be unimportant for conductance quantization in the linear response regime. We 
see that this is a consequence of unitarity. 

Direct backscattering on soft potentials, present in QPC, is usually very small Nevertheless Laughton et al |Tl| l 
suggested a more effective backscattering mechanism: indirect backscattering through localized modes, i.e. standing 
waves trapped in bulges in the channel. They showed numerically that such modes really appear in realistic QPC due 
to overlap of long-range impurity potentials (therefore this process can not be described within Born approximation), 
and that in their presence backscattering is much more intensive, and exact conductance quantization is violated. 

In order to demonstrate this mechanism in our matrix model, regard an extremal case: two nonreflecting scatterers 
(with A = 0) inside the bulge (see Fig.[j]). Bottlenecks arc modelled by two projection operators, P, allowing to pass 
only Np < N transverse modes. We come to expression ([I]) with 
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Tr (f°(E^f°(E)) = Tr j(i - M) * + (PJPt° B ) (l - M) 'll (P a P$) E> f j , (8) 

where M = Ut A (i - P^Utgil - P)i° B . 

If P ^ i interference terms are present, since generally neither f^,Pf^j, nor Pi L is proportional to the unit 
matrix. 

From now then we are concerned only with nonlinearities in I(V s d,V g ) dependence. TV stands here for a number 
of propagating modes. Suppose that the chemical potential of the left reservoir is constant, fiL = E F , while — 

/•./ • ev sd 0- 

Neglecting an energy dependence of scattering matrix at Fermi level, we can then write the dimensionless nonlinear 
contribution to current, i osc = I osc ■ (irh/ 'e) / '(AE F ) as 

i osc (v,u) = J de [n F (e - 1 - v) - n F (e - 1)] cos (eri(e; v,u) - a 2 (e; v, u)) . (9) 

Here A < 1/2 is a numerical constant dependent on details of the system configuration, and the dimensionless variables 
e = E/E F ,ej±{x) = E jX (x)/E F ,v = eV sd /E F ,u = eV g eff /E F ,n = W/{X F /2), and / = L/X F are introduced g§. 

It follows from (^), that the effect can be observed at temperatures of order or less than its characteristic energy 
scale (being an order of magnitude less than the intermode separation, AE±), that is, a few Kelvin. 

In the experiment the rectified voltage was measured, as a response to ac current (FigJ|(a)), 

Vreot(Io,V g ) = l/2% dcblosmcj, R(I smct>,V g ). (10) 
Jo 

This quantity includes effects of both nonlinearity and asymmetry of the QPC response. 
In our approach, though, it is easier to calculate the rectified current: 

Vect(vO) u) = 1/2-k I d(j>v Q sm(j> G(v sm<p, u), (11) 
Jo 

and then define an effective rectified voltage by 

Vr f J t (io,u)/(i R N (u)) = i rect (v a ,u)/(v G N (u)), (12) 

where Gn(u) = R^iu) = N(u)e 2 /irh is the electrical conductance of an ideal QPC with N open modes. 

The results of numerical calculations of v e JJ ct {u) are given in Fig^|. In order to simplify calculations, the contact 
was shaped as a straight channel with impenetrable walls, of (dimensionless) width n. The influence of gate voltage 
was reduced to coordinate-independent lift of the potential instead of changing the channel's width; this agrees with 
experimental situation All the driving voltage was supposed to drop on the right end of the channel. In order to 
evaluate an effect of multiple mode mixing, in the following calculations we have summed up all possible interference 
terms with Gaussian weightes. For example, SG = Yj%=i exp(— (i — j) 2 /2s). Weight factor are introduced 

to reduce contribution of less probable processes with scattering to larger angles (s being a mixing parameter) . 

A qualitative likeness between the experimental and theoretical curves is evident. The period of v e JJ ct (u) as function 
of the effective gate voltage in our calculations corresponds to (^ 0.1 mV), and is thus of the same order as the 
characteristic potential shift in the contact's centre determined in fllCfl ( ~ 100//V). Thus the model gives a proper order- 
of-magnitude description of the effect. The magnitude of observed oscillations (less than 10~ 4 of total conductance) 
can be provided, e.g., by mode mixing on the contact-bank interfaces JHJ or impurity potentials pl]| . 

In conclusion, we have shown that coherent mode mixing in quantum point contacts leads to a new type of 
nonlinearities in this system, which differs from nonlinearity due to conductance quantization and have an energetic 
scale ~ 0.1 mV. These oscillations develop on expence of conductance quantization, and may occur even in systems 
with soft potentials. They are sensitive to the current direction, thus providing an explanation for the observed 
asymmetry of the current-voltage curve in QPC. The magnitude of the oscillations is defined by the backscattering 
rate in QPC governed in the case of soft scattering potentials by indirect backscattering through the localized modes. 
The effect can be observed at temperatures of order AV/ks — IK. 

The effect is a probable cause of recently observed small-scale nonlinearities and asymmetry of I-V curve in QPC. 

The sensitivity of the effect to weak potential changes in the vicinity of the contact suggests that the slow changes 
in time of the impurity configuration near QPC would reveal as a low frequency conductance noise. This may provide 
an explanation for the random telegraph noise observed in the system as well. 
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FIG. 1. Interference effects in a QPC with two nonreflecting scatterers inside a bulge. Localized modes in the bulge appear 
because only Np of N availiable transverse modes can pass the bottlenecks. Inset: Indirect backscattering by a soft potential 
via the localized mode. l)An electron is scattered from the propagating mode a to the localized mode b (standing wave) with a 
small momentum change. 2)It is scattered from the localized mode to the mode c propagating in opposite direction, also with 
a small momentum change. 



FIG. 2. Rectified voltage vs. gate voltage curves, (a) Experimental curve from Taboryski et al. [10]. The solid and broken 
lines correspond to different values of ac current. Inset shows the corresponding conductance vs. gate voltage characteristic. 
The step width is of order 0.1 V, the maximal QPC conductance correspomds to 10 propagating modes, (b) Effective rectified 
voltage vs. effective gate voltage, calculated according to (13), (14). The contact length I = 100/27T, width n — 10.1, vq = 0.03, 
Gaussian mode mixing is supposed with s = 1. The idealized conductance Gn{u) = (2e 2 /h)N(u) is shown as well. The number 
of open modes, N(u), corresponds to the experimental situation, figure (a). 
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